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Abstract 

In this article, we introduce a generalization of the diffusive motion of point-particles 
in a turbulent convective flow with given correlations to a polymer or membrane. In anal- 
ogy to the passive scalar problem we call this the passive polymer or membrane problem. 
We shall focus on the expansion about the marginal limit of velocity-velocity correlations 
which are uncorrected in time and grow with the distance x as \x\'^ , and e small. This rela- 
tion gets modified in the case of polymers and membranes (the marginal advecting flow has 
correlations which are shorter ranged.) The construction is done in three steps: First, we 
reconsider the treatment of the passive scalar problem using the most convenient treatment 
via field theory and renormalization group. We explicitly show why IR-divergences and 
thus the system-size appear in physical observables, which is rather unusual in the context 
of ordinary field- theories, Uke the <j)^-moAe\. We also discuss, why the renormalization 
group can nevertheless be used to sum these divergences and leads to anomalous scaling of 
2n-point correlation functions as e.g. S'^'"'{x) := t) — B(0, t)]^"). In a second step, 

we reformulate the problem in terms of a Langevin equation. This is interesting in its own, 
since it allows for a distinction between single-particle and multi-particle contributions, 
which is not obvious in the Focker-Planck treatment. It also gives an efficient algorithm 
to determine S^" numerically, by measuring the diffusion of particles in a random veloc- 
ity field. In a third and final step, we generalize the Langevin treatment of a particle to 
polymers and membranes, or more generally to an elastic object of inner dimension D with 
< D < 2. These objects can intersect each other. We also analyze what happens when 
self-intersections are no longer allowed. 
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Figure 1.1: Symbolic picture of a turbulent flow. Particles, or equivalently heat is injected in a finite range of size 
L ^ 1/M (dark grey areas), whereas the turbulent flow grows up to scale I ~ 1/m, which finaUy shaU be taken to 
infinity. This is possible, if the total number of particles, or the total heat, injected into the system is conserved. In 
that case, L and not I sets the largest scale in the problem, and multi-point correlation functions with an anomalous 
L-dependence will be observable in a domain of size L, here symboUcaUy shaded in Ught grey. As wiU be shown 
below, they have anomalous corrections depending on L. 

1 Introduction and Outline 

For now more than 5 decades, turbulence has resisted a satisfying theoretical treatment. The 
principle question asked since Kolmogorov's pioneering work [1] in 1941 is whether there 
are corrections to the simple scaling behavior predicted in [1] for higher correlation functions 
[2,3]. The most natural tool to answer this question is the renormalization group. However, 
all attempts to go beyond Kolmogorov's analysis have essentially failed so far. To better pin 
down the problem, simpler toy models have been proposed. The probably most prominent such 
model is the passive scalar model, introduced by Obukhov [4] and Kraichnan [5]. This model 
describes the diffusion of a point-particle in a turbulent flow with given correlations. For sim- 
plicity these correlations are taken to be Gaussian. Nevertheless, the model is far from beeing 
simple, and shows multi-scaling, i.e. higher correlation functions of the particle density scale 
independly of the second moment, characterized by new critical exponents. More explicitly, 
particles, or equivalently heat is injected in a finite range of size L ~ 1/M, whereas the tur- 
bulent flow grows up to a bound of I ~ 1/m, which finally shall be taken to infinity. This is 
possible, if the total number of particles, or the total heat, injected into the system is conserved. 
In that case, L and not I sets the largest scale in the problem, as visualized in figure 1.1. 

In this article, we introduce the generalization from point particles to higher dimensional 
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elastic objects, as e.g. polymers and membranes. In analogy to the passive scalar problem we 
call this the passive polymer or passive membrane problem. 

We start by considering the passive polymer problem. Much has been learned during the 
last years about higher correlation functions due to a common effort of mathematicians and 
physicists [6-27]. Whereas the first to calculate the 4-point function by considering the 0-modes 
of the steady state are [7], the calculatory most convenient scheme, based on the perturbative 
renormalization group, was introduced in [16]. Contrary to the sometimes heard claim, the 
renormalization group is able to handle large eddy motion. The expansion is performed about 
the marginal limit of velocity-velocity correlations which are uncorrelated in time and grow with 
the distance a; as \x\^, and e small, a relation which gets modified for polymers and membranes 
(the marginal advecting flow has correlations which are shorter ranged.) 

The generalization to polymers and membranes is then performed in three steps: First, we 
reconsider the treatment of the passive scalar problem using the most convenient treatment via 
field theory and renormalization group. We explicitly show why IR-divergences and thus the 
system-size appear in physical observables, which is rather unusual in the context of ordinary 
field-theories, like the 0^-model. We also discuss, why the renormalization group can never- 
theless be used to sum these divergences and leads to anomalous scaling of n-point correlation 
functions as e.g. S^"'{x) :— {[Q{x,t) — 0(0, t)]^"). To do so, we determine the full scaling 
dimension of the composite operators 5*^"'™^ := [(V0)^]" [^V©]^™, with \z\ = 1. In a second 
step, we reformulate the problem in terms of a Langevin equation. This is interesting in its own, 
since it allows for a distinction between single-particle and multi-particle contributions, which 
is not obvious in the Focker-Planck treatment. It also gives an efficient algorithm to determine 
S'^'^ numerically, by measuring the diffusion of particles in a random velocity field. In a third 
and final step, we generalize the Langevin treatment of a particle to polymers and membranes, 
or more generally to an elastic object of inner dimension D with < < 2. Our analysis 
will show that the interesting range for e is < £ < 0. For smaller e, the advecting flow 

is irrelevant. For larger e, the polymer or membrane is overstretched. This is the range, where 
already the particle, i.e. the center of mass of the polymer, shows anomalous diffusion. We also 
generalize these considerations to the case of self-avoiding polymers and membranes. 



2 The passive scalar 
2.1 Model 

The advection of a passive scalar field 6 (a; , t) with a; e M'^ the spatial coordinate and t the time, 
is described by the Focker-Planck type equation [4,5] 

[dt + v{x,t)V]Q{x,t) = iyoAQ{x,t) + f{x,t) . (2.1) 

The correlations of the advecting turbulent velocity field v{x,t) are supposed to be Gaussian 
with zero mean and correlations which grow with the distance r as 

{v\x, ty{x\ t')) = Xy^ix -x\t- 1') = Do5(t -t')j^, P'Kk) ,,f^^~lL. > (2.2) 
where 

pi^{k) ■.^6''--— (2.3) 
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is the transversal projector and m some IR-regulator. The dimension of the coupling uq :— 
Dq/uq in units of m is 

e = [uoL ■ (2.4) 
We will see later that e serves as a regulator. Eventually one is interested in the physically 
relevant case of d = 3 and e = 2/3 (Kolmogorov-scaling) or corrections thereto [28]. 
/ is a Gaussian scalar noise with zero mean and correlator 

{fix, t)f{x', t')) = 6{t - t')Gf{Mx, Mx') = 6{t - t')Gf{x - x') . (2.5) 

Often people use Gf{Mx, Mx') = Gf{M\x — x'\), and its Fourier-transform Gf{k/M). How- 
ever we will use the more general case for clearness of derivation. G/ is the source of correlation- 
functions of 9, which otherwise would vanish. Physically, it may be viewed as source and sink 
of tracer-particles, or heat. We will see below when explicitly calculating expectation values 
that Gf sets the largest scale L = 1/M appearing in physical observables. Therefore, we 
demand that G'/(s, s') rapidly decays to zero for s or s' larger than 1. 

The analysis of Eq. (2.1) is most easily done by using a dynamic action [29,30] 

j[e, e, Lv]^ I e{x, t) [dte{x, t) - uoAe{x, t) - v{x, t)ve{x, t) - f{x, t)] . (2.6) 

x,t 

Expectation values are obtained by integrating e""''® ®'^'^' over 6, 9 and averaging over / and 
V. Since v and / are Gaussian, their average can be taken, leading to 

j[e,e]= / Qix,t){dt-iyo^)e{x,t)- f roo^ , i 

Jx,t »'a;,2/,i 

where 

= Q{x,t)Vie{x,t)e{y,t)vMy,-t)K\^-y) (2.8) 
= e{x, v,e(x, t) e{y, v,e(y, t)D'j{x -y)^Bo^ 



2 wvi^ ' 2 



(2.7) 



w\« - -X - -•wv 



e{x,t)eiy,t)Gfix,y) . 

Note that due to the transversal projector in the turbulent interaction, the partial integration from 

e{x, t)Vie{x, t)D'^{x -y)to- ViOix, t)e{x, t)D'^{x - y) is possible since ViT)^{x - 
y) = 0. 

The free response and correlation functions read 

[e{k, uj)e{k', uj')\ = {2n)'^+^S{uj + J)5^{k + k')R{k, u) , R{k, uj) = . _ 



(e(A;, uj)Q{k\ a;'))o = (27r)''+M(u; + uj')5'^{k + k')C{k, uj) , C{k, cu) = 



co^ + {uok^)^ ' 

(2.9) 

where in the last formula Gf{x, x') was supposed to be of the form Gf{x—x'). Most convenient 
is a mixed time and /c-dependent representation 

(e{k, t)e{k', t'))^ = {27r)'^5'^{k + k')R{k, t' -t) , R{k, t) = 0(^)6-'^°'='* . (2.10) 

This also yields the response-function in position space 

R{x,t) ^ e{t) (47rz/oi)"2 e"^ . (2.11) 
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2.2 Perturbative corrections, renormalization of the dynamic action J 

We now study the renormalization of the model, i.e. we want to eliminate all UV-divergent 
terms at £ = 0. It is important to notice that such divergences only come from the insertion 
of the turbulence-interaction ^ 3^ ' but not from the insertion of the source of tracer-particles 
i w. X ^ . To first order in Dq, the only contribution is 



2 



The diagram is without the external legs in the most convenient mixed t and k representation of 
Eq. (2.10) 

[\=h^l^a me-''^^^ (l-^). (2.13) 

In order to clarify the role of the factor 6{t) in . we have to recall that this is an approxi- 
mation for a sharply peaked but nevertheless smooth function around t — 0. This is the reason 
why in Eq. (2.13), we have replaced the 5-distribution by a smoothened one which in the 
limit of ^ will reproduce 5{t). Integrating / dt 0{t)e^^ ^^*'8r^{t) and then taking the limit of 
77 — > thus simply yields a factor of \. Eq. (2.13) becomes 



where Cd is defined as 



{2tiY (A;2 + l)^ 



2r(i + |; 



(2.15) 



r(^)(47r)'^/2 

This leads at leading order to a renormalization of v (denoting with subscript bare quantities) 



where we have introduced a coupling uq and its renormalized counterpart u through 

Do = ZdD (2.17) 
Uq = — = um —- . (2.18) 

We now claim that Eqs. (2.16) and (2.17) are all renormalizations needed, and that even to all 
orders in perturbation theory. Let us first focus on the renormalization of D. There will appear 
diagrams like 

q p q p q p 
> '\A/\» WA#^ > ?^AA»^ -wj* > ^ ^AA«^ > (2 19) 



The passive scalar 



7 



We want to argue that due to the transversal projector in 3^ > this and all similar diagrams are 
finite. Up to an overall factor, and integrating over the time difference between the two vertices, 
they are 



/ 



X 



' + l)][{k + qf] - [q(k + + l){k + q)] 
{k + qY 

+ Q] [jk +pf]-\p{k+p)\ [ip + (fe + P)] 

(A;+p)2 



(2.20) 



Since for large k 



[q{q + IMk + qf] - [q{k + q)][{q + l){k + q)] 
(k + qy 

\p{p + rmk + Pn-[pik + P)][{p + r)ik + p)] 



{k + pf 



0(A;°) 



the integral (2.20) scales for large k as 

d'^k 



I 



(27r) 



d+e 



I 



dk 1 

~kk^ 



(2.21) 



(2.22) 



1. Note that the 
at d = 2 or 4, 



and is thus f/V-convergent for any d and e > 0. This means that 
transversal projectors ensure that no additional divergences appear for e 
since it allows to bring the derivatives always to the external legs. Moreover, no long-range 
interaction can be generated. This argument can be generalized to higher orders in perturbation 
theory, however only the absence of additional divergences for > 2 is immediately apparent. 
We shall not elaborate on this question any longer, since it is not at the center of our analysis. 

Let us now come back to counter-terms for p. By direct inspection, one sees that the only 
diverging diagrams are chains of bubbles, of the form 

^ ^ .... (2.23) 

However, these diagrams are already renormalized by Eq. (2.16). This is easily seen by directly 
summing the perturbative (geometric) series, as e.g. in [31]. The /9-function to all orders in 
perturbation theory thus reads 



u — —eu + I 1 — : — u 
dj 2 



(2.24) 



(Note that we do not take Cd at e = 0; this "minimal subtraction" is completely sufficient at 
1-loop order, but unsufficient for the all order result (2.24).) This /3-function has a fixed point at 

(2.25) 



u 



{d-l)Cd ■ 

One can now define the anomalous dimension 7^ of u as 

dm 



uil- 



1\ Ca 



(2.26) 
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which atu — u* reads (to all orders in e) 

ll^e. (2.27) 

Since G/ is not renormalized, this leads to an anomalous dimension of © (in units of x ~ 1 /m) 
to all orders in perturbation theory as 

77* = -| . (2.28) 

The full dimension of © thus is 

[e>].,f = 1 - 1 • (2.29) 

This simple scaling is only part of the whole story, as we shall see in the next section. 

To summarize: we have constructed a renormalized action which is UV-finite in the limit of 
£ ^ 0, and which gives the IR-scaling for £ > 0. 

2.3 Observables and IR-divergences 

We now want to study correlation-functions as e.g. 

S^^{x-y,t) := ([0(a;,t) -0(y,t)]2«) . (2.30) 

We always choose x and y inside the injection region, thus especially L ^ \x — y\. It will turn 
out that these observables are sensitive to the size L = 1/M of the system, and demand new 
renormalizations. Since from the viewpoint of 0^ -theory this is rather strange, let us study an 
expectation value in the latter theory in order to see where the difference to the passive scalar 
problem lies. Suppose, one wants to evaluate the expectation value 

U{x,y) ■.= ]^{(t>\x)4>\y)) , (2.31) 
for the theory defined by the Hamiltonian in d dimensions 

n[cf>] = I d'x ^{Vcj>{x)f + h(t)\x) . {131) 

(For the difference in definition between Eqs. (2.30) and (2.31) note that for correlation func- 
tions growing with the distance, definition (2.30) has to be used, whereas for decaying correla- 
tion functions, (2.30) is the correct one.) 

Denoting expectation values in the free theory by C(a; — y) := {(f){x)(f){y))Q ~ |x — y\'^~'^, 
the first contributions to U (x, y) are (when setting to self-contractions in the -interaction, 
and neglecting combinatorial factors) 




[/(x,|/)= <>-6<>0*+^' <>0<> + < > +^(^')- (2-33) 



This formulas is to be understood such that the outer points are always x and y and that one 
integrates over the inner points. Since C(s) ~ s^~'^ the term of order h scales as 
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which for for large z becomes 



1 A^z \z\^-^'^ . (2.35) 
It is IR-convergent at least for d close to 4. Now still consider one of the terms of order \? . 

= J d'^s dH Cix - sfC{s - tfC{t - yf . 



(2.36) 



Similar to what has happened in the last section, there is a logarithmic divergence at £ = for 
small s — t, which has to be renormalized. Calculating directly the integral over s — tm the 
regularized theory at d < 4 , this leads to (z :— s — t) 



,0* = jd'zCizf 



-4-d 



A-d 



(2.37) 



where L is an effective IR-cutoff. The question now arises, what L is. Noting that the integral 
over the center of mass (s + t)/2 is IR-convergent with the identical argument that led to Eq. 
(2.35), the effective scale at which the integral (2.37) is cut off, is L = \x — y\. These kind 
of arguments can be continued to higher orders ^ They show three things: First, expectation 
values of physical observables are IR-finite, i.e. boundaries of the system do not enter into 
their calculation. Second, the distances between the observable points set the largest scale L in 
the problem. Third, when varying these distances, L changes and thus the value of diverging 
subdiagrams as (2.37). This gives rise to an anomalous scaling of the observables. The latter is 
most comfortably taken care of by the renormalization group procedure, which also allows for 
a proof of the above statements. 

Let us now turn back to the passive scalar problem, and consider 

S\x - y, t = 0) := ([e(x, 0) - e(y, 0)]^) . (2.38) 
The order contribution is 

/■oo p 

= / dt d'^zd'^z' [R{x-z,t)-R{y-z,t)] 

X [R(x - z', t) - R(y - z', t)] Gf(M\z - z'\) (2.39) 

Using Eq. (2. 1 1) this can be written as 



~ J d<^zdV( [{x-zf+{x-z')^]' [{y-zf+{y-z')^]' " - [{x-z)^+{y-z')^ 



l-d 



r\2 



{y-zy+{x-z') 



l-d 



Gf{M\z-z'\) . (2.40) 



Integrating over both z and z' large, the integral (2.40) scales as 

ds 



J d2'^ss2(i-'i)s-2G'y(Ms) ~ 1 ^ Gf{s) - InL , 



(2.41) 



'a caveat is in order here: Calculating at small but finite values of e := 4 — d, there is always an IR-divergence 
at sufficiently high orders in perturbation theory. To see this take a long chain of bubbles, similar to that of Eq. 
(2.36). The standard way to circumvent this well-known (technical) problem of the massless theory is to take s 
small enough inside each diagram. 
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where the factor of is due to the differences in Eq. (2.40), and L — 1/M is the scale at 
which Gf starts to decay rapidly. S'^{x — y, t) thus explicitly depends on the largest distance 
L > |x — y| in the problem, in the very contrast to the example of the 0^ -theory considered 
above. 

We shall now show that perturbative corrections in D make S"^ depend even stronger on 
L, namely as L"^ at nth order in perturbation theory. To this aim consider the term of order 
D. Using the regularized theory, it appears at two places: First v is renormalized, and thus 
the resulting effective response-function R{x, t) decays faster - this effect would render (2.41) 
IR-convergent. However, there is a second contribution, namely 

D -ciin = j^'^j d'^' / ^ [^(^ - ^) - - ^)] 

X ^ [R{x - z\ t) - R{y - z', t)] Di^ix - y) . {lAT) 

This is only some part of the diagram, and in principle, it has to be closed through a w>. « , 
leading to 

^C^^x . (2.43) 

Since this diagram is again plagued by an IR-divergence, the leading contribution in Eq. (2.42) 
will come from the domain of large z and z' . In that limit, \x — y\h much smaller than both z 
and z' and consequently, 0(x) — 0(y) can be replaced by (a; — y) V6(^i^). Let us again stress 
that this is valid in the domain, where \x\ <^ L and \y\ <^ L, or when using Gf{M\x — y\) 
where \x — y\ <^ L. Due to that replacement, we can now use a very powerful trick: Instead 
of analyzing the IR-divergences of {[Q{x, t) — Q{y, t)]^), or more generally of 5'^"(a; — y,t) :~ 
{[Q{x,t) — G(|/,t)]^") we can analyse the UV-divergences of the composite operator [{x — 
I/) V6(^^)]^'^. The latter however is a standard task in perturbative renormalization. We will 
see in the next section, that this leads to a whole family of operators and anomalous dimensions; 
the operator with the smallest dimension will then give the term which most sensitively depends 
on L. In order to avoid confusions, let us already note that the second moment S'^{r) discussed 
above does not depend on L, since the contribution to the response-function and that in Eq. 
(2.43) just cancel. This can also be obtained exactly [5]. 

2.4 The scaling of S"^"^ and renormalization of composite operators 

As discussed in the last subsection, we now have to study the renormalization of [zV6]^". It 
will turn out that under renormalization this term generates [2;V9]^'^~^^^[(V9)^]. In a sec- 
ond step, [zVQf''-'^z\VQf] generates [zVeY'''-'^z%VQff a.s.o. until also a term of the 
form 2;^"[(V0)^]" is generated. All these operators will mix under renormalization. The eigen 
operator with the smallest dimension will give the term which most strongly depends on L. 
We now treat the general case. Define 

5(n,m) ._ ^2n [(Ve)^]" [^Ve]^'" . (2.44) 

We first observe that the operator product expansion (denoted by ♦ ) is 

D 'wj^ — 



2 ^AA»^ 



(2.45) 
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with 



d 



(ve)= 



n— 1 



+2nmz 



In 



zve)'™ + 2n(n - i)z'"(v*e)(v^e) [(ve)^ 

n-l 



ra-2 



+m(2m - i)z2"z^z^' [(ve)^]'' (zve)^"^-^ . 

Eq. (2.45) then reads 



yo 7 (27rj" fn2_^^2^ — 



(ve)2 - 



(pve)^ 



(2.46) 

+ (^(e") . (2.47) 



Since it!(p, = e ^°P^*©(t), integration over t yields (up to finite terms) 



2UoJ (27r)<^^ p2 ^^2+^2)^ 

T'^ in Eq. (2.46) has the form 



(2.48) 



(2.49) 



Inserting this into Eq. (2.48) and using the formulas from appendix A.l and Eq. (2.15) yields 



u m 

— Grf 

2 e 



1 



d 



1 



d d{d+2) 



(2{BVe){CVQ)+{BC){VQy 



(2.50) 



Specifying T*^ in Eq. (2.49) to its value of Eq. (2.46) gives 



u m 



1 



2 " e d{d + 2) 



[n{d -l){d + 2n + 4m) - 2m(2m - 1)] z'''' [(Ve)^]" [zVO]^"" 
+m{2m -l){d+ l)^2n+2 [(yQ^^n+i [^vej'^^-^j . (2.51) 

The final result when contracting S'^^''^^ with § is 



m 



2 w.^^ 2 £ d(d+2) 



2 £ 



m 



[n((i - l)((i + 2n + 4m) - 2m(2m - 1)] 

m(2m - 1) . (2.52) 



d{d+2) 



This shows that 5(o,n)^ ^(l,n-l)^ ^ ^(n.O) 

mix under renormalization. The first two eigen op- 
erators are 5^"'°^ and 5^""^'^^ := S^^~^'^^ — |»S^"'°), with eigen values 



^(n,o) ^ £ ^ 5(n,o).« C'rf -l){d + 2n) 



2 ^AA»^ — 

D -NAA*] 

2 wx*] 



2 £ 



m 



d(c^ + 2) 



^ ^ ^ ^(n-1,1) 



ii^wCrf l)(d- l)(d + 2n + 2) -2 



2 £ 



m 



d(d + 2) 



(2.53) 
(2.54) 
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More importantly, due to the triangular form of the matrix, the eigen values can just be read 
off from the diagonal. Therefore the eigen operators 5^"'"*^ are multiplicatively renormalizable 
(again denoting by index o bare quantities) 

^(n,m) _ ^(n,m) ^{n,m) 

^(n,m) ^ y _ [^(^ _ + 2n + 4m) - 2m(2m - 1)] . (2.55) 

This yields the anomalous scaling-function ^^^'"^^u) of Sq"''"^^ in units of x as 

dm 
V 1 

= Cd-rr; r Hd -l){d+2n + 4m) - 2m(2m - 1)] . (2.56) 

2 did ~\~ 2 j 

At the IR fixed point u = u* from Eq. (2.25), this is 

^(n,m) _ _ _ ^ ^^^^ _ + 2n + 4m) - 2m(2m - 1)] . (2.57) 

(a — l)(a + 2) 

Taking care of the naive perturbative contribution to the scaling of <So"'"*^ from Eq. (2.28), we 
finally obtain the full scaling-dimension A^"'"*) of S^^'^^ 

^(n,m) ^ ^(n,m) _ 2(n + m)r]* 

= + (n + m)£ 

_ 2n{n - 1 + 2m) m(4m - 1) + + d 

" (dT^) + ' (d + 2)(d-i) ■ ^^-^^^ 

These results have already been obtained in [7,12,16], where only the case m = was given. 
The case m > can be found in [32]. 
The exponents satisfy the inequality 

^{n,m) ^ ^{n-l,m+l) ^ (2.59) 

such that indeed A" := A*^" ") is the smallest of all exponents A*^"^""^'*") and dominates the 
L-dependence of ^^"(r, t) as stated above. Explicitly, 

S'-{r, t) ~ r"(^-) (0^" , A" = -^^^^^^ ■ (2.60) 

Only the second moment does not depend on L. As demonstrated in [5], this is the consequence 
of a conservation law, which allows for an exact calculation of the second moment. 

These results have been tested numerically with different methods in [10,21,23,27]. 

The other question one might ask is why only one of the two factors in Eq. (2.60) depends on 
L, and whether the r-dependence comes out correctly. To understand this point, we recall that 
the first factor is due to the renormalization of u, and thus contributes to the anomalous dimen- 
sion of ©, irrespective of the boundary conditions. The second factor stems from the anomalous 
dimension of the composite operator 5*^"'°^ which was associated to the IR-divergence, i.e. L- 
dependence of ^^"(r, t), and which has two contributions: the proper renormalization of 5'*^" °-' 
as given by Z^"'") or 7("'°), and the renormalization of u; these add up to A("'°) as given in Eq. 
(2.58). Only the combination of these terms contribute to the L-dependence of 5'^'*(r, t). 

Also note that the exponents with m > are also observable, and correspond to observables 
of different symmetries. 
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3 Langevin-description of the passive scalar 

3.1 Model and basic properties 

Let us now turn to a Langevin-description of the passive scalar problem. We start from Eq. (2. 1) 

[dt + v{x, t)V] @{x, t) = uoAe{x, t) + f{x, t) . (3.1) 

Without the source f{x, t), this can easily be converted into a Langevin equation. The question 

arises, how the additional term f{x,t) can be incorporated. We will see below that it corre- 
sponds to the creation and annihilation of particles, and that it can indeed be formulated within 
a Langevin description. However, this is a question of marginal interest, since we had seen in the 
last section that the whole renormalization procedure can be performed without ever specifying 
the correlations of /(x, t), just knowing that they will deliver some IR-cutoff L. We therefore 
start our analysis by studying Eq. (3.1) with f{x, t) = 0. Using standard arguments [33], it is 
transformed into a Langevin equation for the motion of a particle with position r{t) E R"'. Since 
it will turn out later, that to reproduce all expectation values, one has to introduce N particles, 
we will already give the corresponding generalization here. 

dy^it)=v\r^{t),t)-C{t) 
{Ca{t)Cpit')) = 2uo5'^5apS{t - t') . (3.2) 

The vv correlations are the same as in Eq. (2.2). The dynamic action which enforces the 
Langevin-equation to be satisfied reads 

J[r,f,v,C]^Yl [dtfi{t)[dtri{t)-v'{r^{t),t) + Cm] . (3.3) 

a=l •' 

Averaging e"'^''"''''"'^] over C, and v leads to 



J[r,f] = / E [fimri{t) - uofiitf 

■' a=l 



/ 



Do 



AT 



^ a,/3=l 



Mlc c.iWa{t)-rji{t)\ 



f\{t) ■ (3.4) 



Symbolically, this is written as 



„ / AT D \ 

J[r,f] = y dt E {'^oc - i^O'^'^'^a] - ^ E 1 ■ (3.5) 

\a=l Q:,/3=1 / 



Free response and correlation functions are 



C{^) = 5 , C{t) := ^ (Ht) - r{Q)f)^ = uo\t\ . (3.6) 
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3.2 Renormalization of the dynamic action 



In subsection 2.2 we have seen that in the dynamic action (1.1) only vq demands a renomial- 
ization. How does this renormalization show up in the formulation as a Langevin equation? To 
answer this question, write down the first order term in Dq from the expansion of e"-^: 



a,/3=l 



with 



Jk[ra(t)-ri3{t)] 



d+e 
2 



it). 



(3.7) 



(3.8) 



{27vY ' ' (A;2 + m2)^ 

We now have to analyze short-distance divergences, i.e. what happens if the two points come 
close together. This is most easily done using the techniques of multilocal operator product 
expansion, introduced in [34,35], further developed in [36-39] and reviewed in [40]. By the 
dashed line which encircles the two fat points, we indicate points which come close together. 
We now want to express e*^''""^*^"'''^'^*)! through its normal-ordered version :e**^[''"(*)~''/3(*)] :, which 
does not contain any self-contractions. This is 



gifc[ra(t)-r^(t')] ^ .gifc[r„(t)-r^(t') 



e 2d 



(3.9) 



This leads to a drastic simplification: Since ^ (['^a(^) — '^/3(0]^)o equals infinity except for 
q; = /3 for which it vanishes, Eq. (3.8) gives 

d'^k 



where Cd is defined in Eq. (2.15). Symbolically, the diagram is written as 



a/ e 



'm 



d) e 



(3.10) 



(3.11) 



which reminds of Feynman's bra and ket notation. We have added an index m to indicate the 
IR-cutoff. Eq. (3.11) leads to the same renormalization for v as given in Eq. (2.16) for the 
Focker-Planck formulation of the problem. 

Even though equivalent, the treatment in terms of the Langevin equation reveals one impor- 
tant property: The only renormalization of the dynamic action comes from the divergence in 
a single particle trajectory. We will see later, that the renormalization of 5" is due to multi- 
particle diagrams. 



3.3 Simple expectation values and translation table 

In this section, we make more explicit the relation between the two formulations of the problem. 
What we want to calculate are expectation values ofQ{x,t). We first observe that in the limit 

of ^ oo 

1 ^ 

Q{x, t) « BlIx, t):=-Y^ 5\r^{t) - x) . (3.12) 
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If we do not know where the particles started, then obviously 

{e^{x,t)), = 0. (3.13) 

This is also true for higher moments 

{e^{x,tr), = 0. (3.14) 

It is important to note that the limit of — > cxd is necessary in order to suppress correlations 
coming from the same particle, which in the expectation value (3. 14) are of order 
Equivalently in momentum representation, when defining the Fourier-transform as 

Eq. (3.12) will read: 



T^Ji^y'' (3.15) 



1 ^ 

Q{Kt) ^ e^{k,t) = - 5: e-*^^«(*) . (3.16) 

-'^ a=l 

We now check that the free response-function (2.10) of the Focker-Planck formulation is cor- 
rectly reproduced. In contrast to the Focker-Planck formulation, here we have to solve the initial 
time problem explicitly, i.e. suppose that the particles start at position xq at time to. 



Ri^{x-xo,t-to) -.^ {ei^{x,t))J ^ ^ eit-to) 

\ra(to)=xo 



N 

ik[ra{t)-ra (to)] 



E 



e{t - to) 

ra(to)=XO 

Qik[xo-x]Q^^ _ ^^^^ 

ra{to)=xo 



= Rfp{x - xo,t- to) , (3.17) 

as given in Eq. (2. 11). Note that this response-function is a single-particle function, i.e. only the 
response of a single particle « to a change in its (earlier) trajectory contributes. 

Injecting particles with a rate f{x, t) at time t at position x, models the corresponding term 
in Eq. (2. 1), and finally leads to the same correlation function as in the Focker-Planck represen- 
tation 

Ci,{x, t; y, t') = j dto j dt'^ j d'^xo j d'^yo Rl{x - xo,t - to)Ri,{y - yo, t' - t'^) 

X {f{xo,to)f{yo,to)) 

= y dto I d'^.To I d''yoRj^{x-xo,t-to)RL{y-yo:t' -to)Gf{Mxo,Myo) 

= CMx,t;y,t') . (3.18) 

By comparing the terms of order v{x, t) in Eqs. (2.6) and (3.3), we also obtain the equiva- 
lence 

N 

e{x, t)v'e{x, t) ^ e^ix, t)v'eL{x, t) Yl ^ai^y^^Mt) - x) . (3.i9) 

a=l 
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3.4 The scaling of ^S*^" and renormalization of composite operators 

We now have to reproduce the results of section 2.4. This will be done in two steps. First, we 
will show that the OPE of Oi{x, t) := Ql{x, t) V'GlIx, t) = Ea=i r'a{t)6'^{ra{t) - x) with the 
n-th power of ©L(y, t') = J2p=i ^"^{^pit') ~ v) is as in the formulation via a Focker-Planck 
equation encoded in the contraction of t) with a single ©L(y, t')'- 

Ol{x,t)*Qi^{y,tr = nQi^{y,t'r-^ [0\^{x,t) * Qi^{y,t')\ . (3.20) 

The reason is that when contracting fl^{t) in t) with 0L(y, t')'^ singles out a particle a in 

one of the ©l(i/, t')- Correlations of this particle a with rg(t') in another of the factors 6l(|/, t') 
only exist for o; = /3 which is suppressed by a factor of 1 /N. 

Therefore, it is sufficient to show that 0^{x, t) ♦ 0(y, t') is the same in both the Focker- 
Planck and the Langevin formulations. In the Focker-Planck formulation, we have 

0^p(a;, t) ♦ eFp(2/, f) = eFp(a;, t)V'RM^ -y.t-t'). (3.21) 

In the Langevin-formulation, the same expression reads 

TV 1 JV 

t) ♦ eL(i/, t') = E ms^T^it) -x)^-Y. nr^it') - y) 

a=l -'^ 



1 ^ 



JV 




^ifc[ra(t)-a;] ^ Qip[ra{t')-y] 



(3.22) 



where the normal-order sign ":" indicates that contractions between the included operators are 
factored out as e'^^'*^* L This is useful since the normal-ordered operator ■QMrait)-x]Qip[rait )-y] . 
is free of divergences when approaching x and y. For an introduction and review of these 
techniques, see [40]. 

In the next step, the integration variable k is shifted to k ^ k — p 

Partially undoing the normal-order procedure finally leads to 
e^t' -t)j^J2j-^ j-^d^'Py le^'^l^"^*)-^!: :e^^[^«(*')-'^«Wl:e-^''l*-*'le^^[^-^'l . (3.24) 

Note that this is not the standard procedure (as described in [40]) but is particularly useful for 
our purposes. Two routes of argument are now open: The first one consists in the observation 
that in the desired limit of t — > t', :e**'['"«(* )-^«(*)] ; is approximately 1, and thus the integrals over 
k and p factorize, leading to 

e^{x,t)V'Ri^{x-y,t-t'), (3.25) 
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Figure 3.1: The two different types of diagrams, involved in the renormalization of (a) f (left) and (b) 5^" (right). 
In both cases we have drawn the particle trajectories, as well as one interaction. 



which using Eq. (3.17) is the same result as in Eq. (3.21). 

We finally want to argue that the above result becomes exact, when specifying the boundary 
conditions. To that purpose, suppose that we start at time r at position 0. The condition that the 
particle be at position is again expressed as a 5-function and using the Fourier-representation 
^'''{fai'T)) — I (^j^ e*'^"^"^^, we are led to study the expectation value of 

lQilra.iT) f d"'^ ^ik\rJt)-xA _ f _^^^kl\t-T\^-ikx Uilra,{T)+ikTc{t)\ (3 26) 

\ J {2tiY I J {27rY \ /o ■ 

Since ^e*'''"^'^^"'"''^^"^*^^^ = (27r)'^5'^(/c + p), only expectation values of exponentials 



survive for which global "charge neutrality" kj — holds. Therefore, in Eq. (3.24), we 
have to supply an additional factor e"''^''"^^) at initial time r, leading inside the integral to 

.gife[r„(t)-r-„(r)]. .gjp[r„(t')-'-a(t)] • Q-p''\t-t'\-k^\t-T\ _ (3.27) 

(Note that we have already normal-ordered the first exponential, leading to the factor of e"'^^ .) 
The key-observation now is that 

.Qik[ra(t)-ra{T)]. .Qip[ra{t')-ra(t)] . ^ .gife[r-a(t)-ra(T)]gip[r-a(t')-»-a(t)] . Q 28) 

as long as T < i and (by assumption) r < t'. This factorization property is one of the essential 
simplifications for polymers (which are nothing but a random walk), see e.g. section 10 of 
[40]. For our case, it shows that the above stated approximate equivalence between the Focker- 
Planck and Langevin descriptions is indeed exact, as one expects from the equivalence of the 
two equations. 

In conclusion: Since the above arguments show (at least at leading order) the equivalence 
of the both perturbation theories for the renormalization of S'^'^, we obtain the same results as 
in subsection 2.4. 

Let us still give some remarks on the class of diagrams involved in the renormalization of 
5^". Whereas the diagrams which contribute to the renormalization of v are single-particle 
diagrams, i.e. diagrams where one particle interacts with itself, the diagrams which contribute 
to the renormalization of 5'^"(x — y) are multiple-particle diagrams, i.e. diagrams where one 
particle which finally ends at x interacts with another particle, which finally ends at y. 
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L 



Figure 3.2: Intuitive interpretation of — y) via particles being advected by the turbulent flow. Drawn is one 

configuration which contributes to tcy of Eq. (3.38). Note that particles may return into the box, such that one has 
to wait long enough until the probability of return tends to 0. For details cmp. the main text. 



First of all, one can determine the single particle motion, which is super-diffusive. By means of 
a complete RG-analysis, or faster using the method of exact exponent identities [40], we obtain 
due to the non-renormalization of the terms proportional to fdtr and Dq the exact identity 



which is the analogue of Eq. (2.27) ff. 

More interestingly, 5'^" (a; — y, t) can also be obtained in terms of particle trajectories only, 
following [27,41]. We first note that due to Eq. (3. 12), Qh{x., t) is the contribution from particles 
which are created by f{x', t') at time t' <t and position x' and arrive at time t at position x: 



This gives a simple method to evaluate moments of © within a Monte Carlo simulation. The 
following observation helps to render the derivation more transparent: 

For given ^ and v, the process is deterministic. Since C, and v are Gaussian and uncorre- 
lated in time, they are time-reversal invariant. We therefore write Qh{x, t) as integral over all 
trajectories, starting at time i at x and ending at time > i at x', where they are "created" by 



3.5 Interpretation in terms of particle trajectories only 




(3.29) 




(3.30) 



/(x',0: 




a=l 



(3.31) 



where r^{t) satisfies the equation of motion and boundary condition 



dt<{t')^v{rl{t'),t') + Ut), 



rl{t) = X . 



(3.32) 
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We now turn to the evaluation of higher correlation functions. Let us demonstrate the principle 
on the example of the second moment, keeping in mind that finally — > oo: 



[e{x,t)-e{y,t)f 



[ej^{x,t)-ej^{y,t)f 



•I Iv 



1 ^ III r 1 r 



N 



'f'v 



1 ■'^ poo I \ 

— E / dt'((Gf(rS(0,r^(0) + Gf(r^(0,r^(0)-2G'f(rS(0,r^(0))j , 



(3.33) 

where we used that 

(/(r, t)/(r', t')) = b{t - t') GfiMr, Mr') = 5(1 - t') Gfir, r') (3.34) 
is 5-correlated in time. In principle, Gf{x, y) has to fulfill four conditions 
(0 Gf{x,y) = Gf{y,x) . 

(ii) Gf{x,y) — as a consequence of J^. f{r, t) = in order to have no global heating 
(particle conservation). This condition is necessary to reach a steady state. 



(Hi) Gf{x,y) 



. We take it at least exponentially decaying at scale 1. 



(iv) Gf has to be realizable as a stochastic process, such that as a consequence of 



^{f{x,t)-f{y,t))') > one must have G f{x,x) + G f{y,y) - 2G^x,y) > 

A possible choice for Gf{x, y) that satisfies the above conditions is 

Gf{x,y) A,A,e-^'-^' = 4e-^'-^'(2x^ - d){2y'' - d) . 
Constraint (iv) is satisfied since 

Gi{x, x) + Gf{y, y)-2 Gf{x, y) = A Q-K'^'+y') [d (f' - e^') + 2 e^'x^ - 2 e^'y^ 



(3.35) 



> . 
(3.36) 

In simulations, people [27,41] use Gf{x,y) = &{\x — y\ < 1) (with being the step 
function), which violates condition (//), but appearantly leads to good results. (Particle con- 
servation seems to be no problem, when working at fixed particle number. For another argu- 
ment see [27].) As explained in section 2.1, we make the ansatz that Gf{x,y) depends both 
on X and y and not only on the difference x — y. Similar in spirit to [27,41] would thus be 
Gf{x,y) = &{\x\ < l)&{\y\ < 1). Eq. (3.33) then acquires the form 

[e{x,t)-e{y,t)f 



N 



7V2 



a,f3=l 



E / dt' {( &iK{t')\<L)-e{Kif)\<L) 0(|r^(t')|<i.)-e(|r|(t')|<i^) 



(3.37) 
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We see that this is the I'-average of a quantity which can be interpreted as a time. Let us in 
general define 



a,/3=l * 

(3.38) 

Then (jO(x, t) — Q{y, t)f^ is the v-average of txy'. 



p{x,t)-Q{y,t)Y) = {t,yl . (3.39) 
Analogously, for the fourth moment we can write using that / is Gaussian and — > oo 

[e{x,t)-e{y,t)t)^-- j dt' dt" 

GfK{t'),r^^{t')) + Gf{rl{t'),r${t')) - 2Gf {rl{t%rl{t')) 
Gfir^if'Ulit")) + Gf{ry{t"lrl{t")) - 2Gf {r^^{t"),rl{t")) 



X 



^3{{t.yf \ . (3.40) 

In general, the 2n-th moment is 

[Q{x,t)-QMf-)=^^{{t^yTl . (3.41) 

One can now easily check the deviation from Gaussian behaviour by analyzing connected ex- 
pectation values. Eq. (3.41) also tells us that the anomalous behaviour of S'^" comes from rare 
large events, i.e. rare large t,j.y. Moreover, the measured ensemble of tj-y not only contains the 
information about the moments, but even the complete probability distribution function. 

Let us give some comments on simulations [41,27]: In the above formulation, we have 
always averaged over replicas for the particles and their corresponding thermal noises, finally 
taking N ^ oo. For the 2n-th moment, it is indeed sufficient to keep 2n different particles. 
Moreover, the thermal noise can be dropped for particles which are not starting at the same 
point. Finally, by going to the relative coordinate system, and using G/ of the form Gf{x — y), 
only the 2n — 1 relative coordinates have to be propagated. 

Also note that with a little bit of work, using Eq. (3.29) and the anomalous dimension of 
(V©)^", one can again obtain the anomalous scaling behaviour as given in Eq. (2.60). 



4 Generalization to polymers and membranes 

4.1 Construction of the generalized model 

We are now in a position to generalize the above considerations to polymers and polymerized 
(tethered) membranes. 

To this aim, we introduce a polymer or polymerized tethered membrane [42,43] with coor- 
dinates 

X e — > r(x) e , (4.1) 
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where we think of D between and 2, and in particular of D = for a particle, D = 1 for a 
polymer and D = 2 for a membrane. For polymers, x measures the length along the polymer; 
for membranes, x belongs to a 2-dimensional coordinate system. r(x) is the position of the 
monomer x in imbedding space. The standard model for polymers is due to Edwards [44] and 
reads generalized to membranes [42,43] 

7i[r] = j d^x^(Vr(x))2 + bjd^^j d^y 6\r{^) - r(y)) . (4.2) 

The second term punishes self-intersections of the membrane, making the membrane self- 
avoiding. In what follows, we shall study phantom membranes, i.e. drop the term proportional 
to b. We shall discuss the general case in section 4.5. 

As for particles, we introduce A'^ copies, labeled by lower Greek indices. For the a-th 
polymer, the equation of motion for the monomer x at time t and with coordinates ra(x, t) then 
reads 

atr^(x, t) = A.r^(x, t) + ^^(r«(x, t),t) + C(x, t) . (4.3) 

Note that the elasticity of the polymer or membrane has been scaled to 1. The w-correlations 
are the same as in Eq. (2.2), and we already note that e = is not the marginal case. The 
thermal noise is Gaussian with zero mean and correlations 



C(x, i)C^(y, t')) = 2uoS''5^05it - t')5^(x - y) . 



(4.4) 



The dynamic action which enforces the Langevin-equation to be satisfied reads in analogy to 
Eq. (3.3) 



N 



a=l 



J[r,r,v,C] = E / dtd^xf=;(x,t) air;(x, t) - A.r;(x, t) - t;^(r„(x, t), t) - C(x, t) 



Averaging e •^I''''''"'^] over ( and v leads to 



(4.5) 



J[r, r] 



N 

dtd^x E \f-ai^,mrai^,t) - f^(x,t)A.r^(x,t) - m-ai^,t)' 

a=l 

gife[ra(x,t)-r-^(y,t)] 

(A;^ -|- vn?)^ 



N 



■^/dtd-xd^y E r],M 

a,/3=l 



/ 



d'^k 

(2^ 



P'^{k) 



fUy.t). (4.6) 



Symbolically, this is written as 



J[r, r] 



]]-^jdtd-^d-y E „^ 



where we used the abbreviations 



(4.7) 



= ^a(x, i)9tr«(x, t) , = ^a(x, i)(-Ax)r«(x, t) , = f„(x, tf 

AdU „jfe[r-a(x,t)-r-;3(y,t)] 



rvw^ = r^(x,t) 



{k^ -\- rn?) ''2^ 



(4.8) 
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The free single-membrane response- and correlation-functions are 

1 



C{k,uj) = 



, i?(k, t) = e(t)e-*'^ , i?(x, t) = e(i)(47ri)-^/2g-xV4t 



(4.9) 



At equal times, the free correlation-function has the simple form 

c(x-y,o):=l(Hx,<)-.(y.*n=p3^ 

All other free correlations vanish. 



x-y 



2-D 



(4.10) 



4.2 Renormalization of the dynamic action 

As in section 3.2, we now have to analyze possible renormalizations of the dynamic action. As 
in Eq. (3.7), we start from 

n ^ 



a,/3=l 



(4.11) 



which is more explicitly 



d^yr;(x,t) 



Ad I. pife[ra(x,t)-ra(y,t)] 



f^(y,t). (4.12) 



We will see below that contracting a response-field will give no contribution. We are thus left 
to normal-order the r.h.s. of Eq. (4.12) by virtue of a generalization of Eq. (3.9) 



gifc[r„(x,t)-r^(y/)] ^ .gife[ra(x,t)-r^(y/)] . g-fc'^([ra(x,t)-r^(y/)]2)^ 

Only the term with a — (3 gives a non-zero contribution: 



(4.13) 



/d^yf;(x,t) 



(2^ 



.gifc[ra(x,t)-ra(y,t)] . (2-D)S£, '^^l^-yl 

(/c^ -|- m?)^ 



ri{y,t). (4.14) 



The leading UV-divergence is obtained by expanding : e**^[''«(^'*)~''«(y'*)l : and keeping the term 
of order 1 . (Note that such an expansion is justified since the normal-ordered product itself does 
not contain any divergence.) Using the rotational invariance of the integration over k yields 



(4.15) 



Performing the integrals over y and k finally gives the contraction of L J and projection onto 



'm 



D 



m 



D 
' 2-D 



(4.16) 
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2D 

5 = e + — (4.17) 

- jitd 1(2 - D)So\ (4^)./.r(fe-)r(.) 

t=XpL.y(2- D)S„\^ . (4.18) 

This is the generalization of Eq. (3.1 1) to membranes. We have put the index m to the diagram 
to remind the reader that m acts as regulator. Note that now the dimensional regularization 
parameter is 5 instead of e. 

Other renormalizations for the dynamic action do not appear: First, contracting in Eq. (4.12) 
one of the response-fields, say r^(x, t) leads to a factor of /c\ which together with P'^^{k) gives 
0. This argument is generalized to all orders in perturbation theory upon remarking that the same 
factor always appears in the MOPE for the interaction which is most advanced in time. One 
can also show that the emerging diagrams themselves also vanish, due to a response-function at 
equal times. 

Let us now turn back to the renormalization-group functions. Bare and renormalized cou- 
plings are 

u^^^^u^^m' (4.19) 

,,'i-D /72-D 

Z. = l-'^il--^^. Z^ = l. (4.20) 



d 

In analogy to Eq. (2.24), the /3-function reads 



u = -6u + f 1 - 3 ) + 0{u^) . (4.21) 



2-D \ d 



It has a non-trivial IR-attractive fixed point at 



U- = ^Z3±S. (4.22) 



We can as in Eq. (2.26) define the anomalous dimension 7,^ of u as 

7^(1*) := m^lnZ^ 
dm 

= .(1-5)^. (4.23) 

which atu — u* reads 

it = . (4.24) 
A new exponent is associated to the equal-time inner-membrane correlation function 

C„.(x-y) :=^((r(x,t)-r(y,t))^ 

\2k 



^ ■ (4.25) 
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From Eq. (4.24) we obtain the result to order 5 



2 , J. (4.26) 

Using that neither the term / fr, nor / f(— A)r, nor the interaction proportional to Dq = D is 
renormalized, the scaling dimension of r can as in Eq. (3.29) be obtained exactly: 



((r„(x, t) -ra(x,t'))^ 



Since ~ i 



[ra{x,t) -ra{y,t)y 



~ X 



|2fc* 



\t-t'\^- 



K — 



(4.27) 



(4.28) 



2-D 



In the limit of 5 — > 0, this reproduces Eq. (4.26). It is important to note that for £ > 0, 
the exponent k* is larger than 1, thus the membrane over-stretched, and the description of a 
membrane via Eq. (4.2) with only a harmonic elastic term 7i[r\ = J d^x i[Vr(x)]^ breaks 
down. This coincides with the range of s, for which already a single particle, hence the center 
of mass of the membrane, exhibits anomalous diffusion. This range will not be described by 
our model. In experiments one has indeed observed destruction of polymers by a turbulent flow 
[45,46]. 



4.3 Higher moments: the scaling of .S*^" 

Let us now address the question of higher correlation functions. To this aim, we first have to 
generalize the particle-density Qi,{x, t) defined in Eq. (3.12) to polymers and membranes. Be 
V :— J d^x the volume of a membrane, then define 



1 ^ r 



(4.29) 



Note that Qmix, t) has a well-defined limit both for ^ oo and V — > cxd. Even though we 
average over all monomers, it will later turn out that Qniix, t) can be interpreted in terms of a 
single monomer, in analogy of the discussion in section 3.5. 

Let us again consider S'^''''{x, y) := ^[0m(ic, t) — Qm{y, ^)]^")- As in Eq. (2.45), we have to 
study the contraction ofT'^V'QmiV: t)V^0m(y, t) with 



Jint ^/ dtd^xd^x' \j-^,P"'ik) 



-,ik[ra{x,t)-ris (x',t)] 



This reads 

1 f I^D^ f,n^^ ( ^\p-k) /d^([±AO 



^ a=13=l^ ■' 



NV 

xT\p-k)\i + ky 



(27r)' 



(27r)^ 



7,5=1 



♦ ^ / drd^xd^x' E r-(x,i-r) [/^i^'^^(^) 



J,k{r^ {x.,t-T)-rs (x' ,t-T) 

(k^ + rn?)^ 



(4.30) 
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where we have already shifted the I- and p-integration for later convenience. Following the 
partial normal ordering procedure of section 3.4 and noting z = x — y, z' = x' — y' leads to 



xT'\p-ky{k + iy 



{p-kYP''\k){k + lf 



X ■Qi{v~k)[ra(y,t)-ra{yi,t-T)\. ^ip[ra{x,t-T)-x] ^i{x-y){p-k)-(p-k)^C{z,T) j^^^ ^■j 

^ .gj(fc+0[r-/3(y',t)-r^(x',t-T)] . ^il[ri3(y,t-r)-x] ^i(x-y)ik+l)-{k+iy^C(z' ,t) j^)^^/ ^ _ (4.31) 

This expression allows for one exact simplification, namely 

(p - kyP^^k) {k + If = p''P''\k) I'' . (4.32) 
Similar to what has been done after Eq. (3.24), we use the approximations 

.gi(p-fe)[ra(y,t)-ra(x,t-T)] . ^ 

.gi(fe+0[r;3(y',t)-r;3(x',t-T)] . ^ ^ (4.33) 

which were exact for the particle (see the discussion in section 3.4). 

Next, since we are searching for the leading pole in S = e + which comes from the 
region, where k becomes large and where r, z and z' become small simultaneously, we can 
replace k + lhy k and p — khy —k. This leads to 

r^v*e„,(?/,t)v^0„,(?/,t)* Jint 

/ 1 s2 iV JV „ jd ]di ]dl. n"" P"-^ ( h\l'^ 



d+e 



^-Bo(^) EE /d^ydVd^zd^z'dr^y^^T'W^^ 

^gip[ra(x,t-r)-a;] ^il[rfi{y,t-T)-x] Q-k'^[C(z,T)+C{z' ,t)] j^^^^ t)R{z', t) 

= Do /d^z /d^z' [dr [ ^ T-k^k^ /"'^^U e-^^[^(---)+^(-^-)li?(z, r)R{z', r) 

X V"e„,(x, t - r) V^Oraix, t-T) . (4.34) 

One now has to evaluate the integral over k, z, z' and r, with m as an IR-regulator. The 
calculation of the leading pole in 5 is substantially simplified by moving the regulator from the 
/c-integration to the r-integration. For the leading pole, we have (setting vq — 1 for calculational 
convenience) 



-4 



Do/d^z/d^z'^" " dr j^J'^k^yt_^e-''^^^^^^^^ 

(4.35) 
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where we have abbreviated := V"6m(a;, t — t). Since the integral scales like m ^, it can 
equivalently be written as 



5 dm 
_ 4Do m-^ 
~ 2-D 5 



-4 
.2-D 



dr 



d u TDab 



i?(z,r)i?(z',r)e^e^ 



(4.36) 



In order to proceed, we specify T"^ as 

The integral over k can still be performed, leading to 

4 m^^ r(l-l) 



(4.37) 



Dn 



2-D 6 r(f)(47r)'^/2 



^— /d^z /d^z' [C(z,l) + C(z',l)]^-^i?(z,l)i?(z',l) 



X 



^^(vej^ (i-l) + -^(Bc)(ve„)^ - -^(ve„B)(ve„c) 



d{d+2) 

We now introduce the abbreviation 

1 8 r(i-f 



(4.38) 



d^z / d^z' [C(z, 1) + C(z', 1)] i?(z, l)i?(z', 1) , (4.39) 



■ Cf 2-L'r(f)(47r)'^/2 
which is understood to be evaluated at vq = 1. In the limit of 5 it reads 

4r( 



tD 



,2-D 



D 
2-D 



d^z j d^z' [C(z, 1) + C(z', 1)]"^ i?(z, l)i?(z', 1) 
dz' 



Sd[{2-D)Sd\-^t(,^ 

^^^"^ ^l-z'' [ ^{^r [C{z, 1) + C{z', 1)]-^ i?(z, l)R{z', 1) 

ro J Z J Z 



(2 - D) [(2 - D)5d] 



(4.40) 



Note that only depends on D, but not on d. For D — 0, the integrals over z and z' get 
localized at z = z' = 0, and since in that limit C(z, 1) — 1 and R{z, 1) ^ 1, 1 and we 

recover our earlier results of section 2.4. For general D, is smaller than 1. Its explicit value 
is plotted in figure 4. 1 . 

With this abbreviation, Eq. (4.38) becomes 



cyi 



Do m-^ 



Setting 



D tD 



d+l 



dj ' d{d+2)^ ' ^"'^ d{d+2) 



(ve„,s)(ve^c) 



^{n,m) ._ ^2n 



(4.41) 
(4.42) 
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Figure 4.1: as defined in Eq. (4.40). The crosses have been obtained numerically, the soUd line interpolates 
between theses values. 



yields for T'^ as in Eq. (2.46) 



d 



ji— 1 



2n 



(ve. 



2m 



n-2 



\2m 



.,2n 



^'v^e^n + v^e^n (ve 



n—l 



2m- 1 



2m-2 



-\-2nmz 

Insertion into Eq. (4.41) gives 

Do -Ci^/^— ^ { [n{d - l)(d + 2n + 4m) - 2m(2m - 1)] 5^^) 



(4.43) 



2 5 



+m(2m - l)(d+ . 
Reestablishing the necessary factors of z/q gives the final result 



(4.44) 



m 



1 



2""" (5 (i((i+2) 



_|_ ^(n+l,m-l)^^DjD"^ 



- + 2n + 4m) - 2m(2m - 1)] 

-r ^) 
^ d+l 



6 d{d+2) 



m(2m — 1) . 



(4.45) 



Eq. (4.45) is formally equivalent to Eq. (2.52) upon replacing Cd by C^I . As in Eq. (2.55) 
this yields for the eigen-operators S^ 



(n,m) 
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Zin,m) _ ^ _ ^ ^^^^ ^ [n{d -l){d + 2n + 4m) - 2m(2m - 1)] (4.46) 
^ did ~\~ 2ij 

and with the help of Eq. (4.45), we can evaluate the anomalous exponents 7('*'"*) as defined in 
Eqs. (2.56) and (2.57): 

d 



2~"" d{d+2) 
At the IR fixed point u — u* from Eq. (4.22), this reads 

5/^ 



n{d - l){d + 2n + Am) - 2m(2m - 1)] . (4.47) 



/m 



(c^- l)(d + 2) 
The full dimension A^'"*) of the operator 5*^" "*^ then is 
Aj^'"*) = (n + m)£ + T^,"'"*) 
= — (n + mj 



- + 2n + 4m) - 2m(2m - 1)] 



(4.48) 



2-D 



-\-5 \n-\-m 



61 



D 



{d-l){d + 2) 



[n{d -l){d + 2n + 4m) - 2m(2m - 1)] (4.49) 



The term {n + m)e is equivalent to the corresponding term in Eq. (2.58). It will be derived in 
the next section. The contribution to 5"^" is due to the term for m = 



^(n,0) 



n 



2-D ' ' \~ ~ d + 2 ^ 
This gives the final result for S^'^^\r) with r := \x — y\ 

r\-[f^+5(i-^"^)] 



(4.50) 



(4.51) 



Note that for D > 0, < 1 and already the second moment (n — 1) has an anomalous 
contribution at order 5 

A(^'°) = +5(1-1^). (4.52) 



4.4 Physical interpretation 

In section 3.5 we have interpreted equal-time correlation-functions of moments of 9 as ex- 
pectation values of moments of the time t^y which is constructed from the motion of particles 
starting at x and y. This discussion in terms of particle-trajectories was quite general and can 
immediately be carried over to polymers or membranes. In generalisation of Eqs. (3.31) and 
(4.29) we write 



1 r foo 



(4.53) 
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Again, it is to be understood that (x, t') satisfies the equation of motion (4.3) and that (x, t) 

I a 

1 1 ^ 



X. Then define t^y as 



a,/3=l ■ 

(rS(z, t'),r^^{z', t')) + Gf (r^(z, t'), r|(z', t')) - 2Gf (rS(z, t'), r|(z', t')))^ 

(4.54) 

Note that the expectation value is independent of z and z' such that the average over z and z' 
can be dropped. We thus can alternatively define 

(Gf (rS(z, 0, r^(z', i')) + (r^(z, i'), r|(z', i')) - 2Gf{rl{z, i'), r|(z', i')))^ • 

(4.55) 

This object is the analog of Eq. (3.38), with a single particle a replaced by an arbitrarily chosen 
mononmor z on the membrane a. 

Let us now come to the evaluation of observables. In Eq. (4.27), we have seen that the 
scaling of time and space is related by t ~ r^~^, such that 

{e{x, t)e{y, t)) r^\x- y\^-' =\x- y\^-^ . (4.56) 

This result is identical to the particle case, such that without any proper renormalization, ( V0)^" 
would have dimension —ns, establishing the first term in Eq. (4.49). Eq. (4.56) also yields the 
non L-dependent term in Eq. (4.51). 

Finally, in a computer-experiment, one can then measure (r := \x — y\) 

5(-)(r) = {mx,t) - eiy,tr) A^-^) . 

(4.57) 

Let us state explicitly the result for polymers, which are probably easier to simulate than mem- 
branes. Using that the integral in Eq. (4.40) for = 1 gives = 0.456143, we obtain 

, , , , X ji[-2+5('l-0.456143^^±§|i)l 

4S.e.(r) = ([e(x, t) - e(,, tr) ^ r"(-^) (0 ^ ^ ^ {{t^yTl . 

(4.58) 

We also want to give some practical hints: 

One should best use a box of size L and let particles or monomers start at position ±r/2, 
with r <^ L, where the box extends to L/2 in every direction. Also note that the evaluation of 
the 2n-th moment demands to propagate 2n membranes, a number which may be reduced to 
2n — 1 when going to the relative coordinate system. 

Another point is whether one should best measure the L- or the r-dependence. For particles, 
the best observable is the L-dependence, since it directly gives the multi-scaling exponent. For 
polymers, one may instead measure what happens when r — > Ar and L — > X~^L, since for that 
rescaling 

- . (4.59) 
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4.5 Self-avoidance 

Physical membranes are always self-avoiding, i.e. they are described by model (4.2) with 6 > 0. 
One would therefore like to have a combined treatment of self-avoidance and passive advection. 
In general, such a treatment is impossible, since the upper critical dimensions are different and 
not for both couplings exists a small control parameter. In the present problem however, the 
control parameters, i.e. dimensions of the couplings are 

^ — mm = JZTd + ^ 

7 := [bo]^ ^2D- , (4.60) 

where we have introduced a scale ~ 1/x which is common for self-avoidance and which is 
related to m by 

2-D 

m — fj,2 . (4.61) 

By choosing the range e of the turbulent advection and the dimension of imbedding space d, 
both S and 7 can be set to zero. This is the common expansion point. Note that this expansion 
point is for 6 = — |, such that the advecting turbulent field is indeed long-range correlated. 
Short-range correlated turbulent disorder is principally different, since under renormalization it 
generates potential disorder, and physics is described by a new universal fixed point. At least 
this has been observed for static disorder in [39]. 

Let us now turn to a diagrammatic analysis. The most complicated diagrams come form 
the correction to self-avoidance by the turbulent advection and vice versa. First of all, the 
turbulent advection is long-range correlated and thus not corrected by self-avoidance, which is 
short-range. On the other hand, the turbulent advection can correct self-avoidance through the 
diagram 

/ -VVA.^ f9\ 

\ wA» \ ♦yw 

where we have denoted with 



(4.62) 



^ . = 2 j d'^^j d^y j ^^f(x,t)(zA;)e*'^l'-("'*)-'"(y'*)l (4.63) 

the self-avoidance interaction. Explicitly calculating the diagram shows that due to the transver- 
sal projector in the turbulent advection, the MOPE coefficient 

.) (4.64) 



v\A* -•^AAJ 



identically vanishes. This substantially simplifies the analysis, since now only the correction 
to the fields and to the scaling of time intervenes in the RG-analysis. We therefore introduce 
renormalized fields r, f and a renormaUzed time. To this aim, we replace throughout this article 

r — > ro and f fo, as well as t — > X^t. This procedure is more formal than that employed 
in the rest of this article, but necessary to obtain the renormalization group /3-functions. Define 
(again denoting with subscript bare quantities) 

To \fZr 

fo = VZf 

Ao = Z,,\ . (4.65) 
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This gives the dynamic action (with summation over the replicas and integration variables sup- 
pressed, setting also uq — > 1), for which we first give the bare, and second the renormalized 
version: 



J[r, r] 



•wv. ~t~ ^0 '^0 '^^^^^^ 



ZZ wv*- -)- \l Z Z Z\\ ^AA;^ — Z\Z\ wv^wv — Zn "wx* •vw -)- bXZf) wv* 



(4.66) 

Since there is no counter-term for , 

ZZ^l. (4.67) 
Second, from the term proportional to -4- which is corrected by self-avoidance, we obtain 

Z^ = Z^^ 



i--{2-d: 



(4.68) 



where we have used the static notation [36,40], since the dynamic diagrams involved in the 

renormalization of self-avoidance can all be reduced to static ones [37,40]. We recall the nota- 
tion (see e.g. [40]) that (|Q -t-^ means the residue of the pole in I/7 of the diagram -t^^^. 

The latter is defined as the integral over the MOPE-coefficient |-t-^ , cut off at scale L. The 
diagram and Z-factor is as defined in [36], where a different normalization was used. However 
the final result is only sensitive to the ratio of diagrams; moreover since as discussed above 
there is no diagram mixing b and u, the result only depends on the ratio of the two diagrams 
involved in the renormalization of pure self-avoidance. We shall therefore in the following treat 
overall normalizations rather sloppily. 

The term proportional to is only renormalized by the turbulent advection [37,40] and 
reads 

1m/ 



Z^Z 



Z \'vy>>i #wv 



These relations can be solved for Z 



Z 



ZZ. 



l--(2-D)(i 



15 



(4.69) 



(4.70) 



The relation between bare and renormalized couplings for the turbulent advection are in the 
absence of a proper renormalization of u, i.e. Zu — 1 



Uo 



um'Z-'Z-^/'Z^' 



um 



um'z'-^/'Z^' 



1 - -D 

7 



^ + 



u 



5 2, — D \'^^^ 



(4.71) 



leading to the /9„-function 

d 



dm 



u = u 



-S-b- 



2D 
'2-D 



1 

2 — D X"^/^ ♦y^ 



^AA♦V^A, ) 



(4.72) 
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Figure 4.2: Flow of Eqs. (4.78) for J = 1, and 7 — 0.6, 7 = 1.5, 7 = 2.5 respectively. 



Self-avoidance is renormalized by 



^6 = 1 + - 

7 



6 /^-, 



60 and h are thus related by 

h 



= 1 + 



7 



This gives the /Sb-function 

d 



b2-D 
772 



d 



(3b{b,u) := // 



=6 



-7+6 (i 



■7 2 



(4.73) 



M(i + 2 / 



(4.74) 



/ 2 4 XvyS^^wv 



(4.75) 



Eqs. (4.75) and (4.72) determine the critical point, and (4.70) then gives the size-exponent 

InZ (4.76) 



(I. \ 2-D Id 



2-D 



1-b 



(4.77) 



at this critical point. 

For the remainder of this section, and to allow for a simpler analysis, we will specify to 
polymers. We also use our freedom in reparametrization of u to set 
gives the system of equations 



w*v>/v \ = i_ This 



(3u{b, u) — u [—S + b + u\ 
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b -7 + 26+ -M 



(4.78) 



K,{b, u) 




We can distinguish three fixed points: 



The pure turbulence fixed point u* 



5 and 6* = is stable for 7 < 



\5. The value of k is 



(//) The pure self-avoidance fixed point b* — ^ and = is stable for 7 > 25. The value of 



Note that (/) and (//) reproduce the result of the preceding sections and for self-avoiding poly- 
mers respectively. This completes the discussion of properties of a single self -avoiding polymer 
(or membrane) in a turbulent flow. 

The next question is how multiple membrane properties, especially the scaling functions 
5"^" are modified. Two routes may be taken: either one considers real physical membranes 
which are mutually self-avoiding. However, then already the expectation value of Gm(a;, t)^ 
would vanish, since never two monomers can arrive at the same point x at time t, due to self- 
avoidance. Interesting expectation values are ©^(a;, t)0i„(y, i)- In the case of no turbulent 
advection, they are known as contact exponets [47,35,40]; they are also related to the scaling 
dimension of operators in scalar field-theory. 

The other possible generalization, which is less physical, is to impose self-avoidance only 
between monomers of the same membrane. Then, the diagrams evaluated in section 4.3 are 
complete and one simply has to use the modified expressions for k* and for the fixed-point 
value of u in order to obtain the new multiscaling exponents. 

It should also be possible to take into account the back-reaction of the membrane on itself, 
at least approximately. A treatment a la Zimm [48] would in generalization to [38] lead to a 
triple ^-expansion, namely ind — A, D — 1 and 5. 

We leave the exploration of these ideas for future research. 



In this article, we have shown how the multiscaling found in the passive scalar problem car- 
ries over to extended elastic objects as polymers and membranes. This was possible by first 
reformulating the problem in terms of the advection of particles, which then allowed for a gen- 
eralization to polymers and polymerized membranes. Similar to the passive scalar case, we 
have calculated the anomalous exponents to first order in a perturbative expansion. We have 
also discussed, how these quantities can be measured numerically, by studying the drift of two 
monomers, which sit on different polymers or membranes. 



(///) 



The mixed fixed point b* = I7 — |5, u* 
of K at the fixed point is = | + ^ + 



|(5 - I7 is stable for |(5 < 7 < 2S. The value 



J- 

10- 
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A Appendices 
A.1 Some integrals 

For a (i-dimensional rotationally invariant integral over p, which may be = / d'^p e~ ^<=i , 
we have with some constant C 



In the case that all Aj equal A, the latter reads with the same constant C 



j d'^pe-^p' ^CA-'^/^ . 



Using this, we obtain by differentiating with respect to Aj moments of Pi, e.g. 



1 

d 



{p 



2\2 



n -1 



2 2 
P1P2 



d{d + 2) 
1 



d{d + 2) 



Another intelligent way of doing this is to study 

Jp 

This yields e.g. 



(A.1) 



(A.2) 



(A.3) 
(A.4) 
(A.5) 



(A.6) 



jy?\ j{ap){hp){cp){dp) = ^^^^ [{ah){cd) + {ac){hd) + (a(i)(6c)] , (A.7) 

where the global prefactor is most easily checked by setting a = h = c = d and comparing with 
Eq. (A.4). We need 



{xp){yp){lpf = ^^^^ {{xy)l^ + 2(/a;)(/|/)) . (A.8) 
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